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Abstract. 

Motivated by the possible conflict between the Navarro- Frenk- White (NFW) model 
predictions for the dark matter contents of galactic systems and its correlation with 
baryonic surface density, we will explore an alternative paradigm for the description 
of dark matter halos. Such an alternative emerges from Tsallis' non-extensive 
thermodynamics applied to self-gravitating systems and leads to the so-called "stellar 
polytrope" (SP) model. We consider that this could be a better approach to real 
structures rather than the isothermal model, given the fact that the first one takes 
into account the non-extensivity of energy and entropy present in these type of systems 
characterized by long-range interactions. We compare a halo based on the Navarro- 
Frenk- White (NFW) and one which follows the SP description. Analyzing the dark 
matter contents estimated by means of global physical parameters of galactic disks, 
obtained from a sample of actual galaxies, with the ones of the unobserved dark matter 
halos, we conclude that the SP model is favored over the NFW model in such a 
comparison. 
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1. Introduction 



The standard statistical mechanical treatment for self-gravitational systems is provided 
by the micro-canonical ensemble, leading to a regime dominated by gravity and 
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characterized by well known effects such as the non-extensive nature of energy and 
entropy, negative heat capacities and the so-called gravothermal instability. An 
alternative formalism that allows also non-extensive forms for entropy and energy under 
simplified assumptions, has been developed by Tsallis ] 1\ and applied to self-gravitating 
systems [21 El E] under the assumption of a kinetic theory treatment and a mean field 
approximation. As opposed to the Maxwell-Boltzmann distribution that follows as the 
equilibrium state associated with the usual Boltzmann-Gibbs entropy functional, the 
Tsallis' functional yields as equilibrium state the "stellar polytrope" , characterized by a 
polytropic equation of state with index n (see equation (J2J)). The stellar polytrope (SP) 
yields a Maxwell-Boltzmann distribution function (the isothermal sphere) in the limit 
n —>■ oo. This index is related to the "non-extensivity" parameter q of Tsallis entropy 
functional, so that the "extensivity" limit q — > 1 corresponds to the isothermal sphere. 

Although the self-gravitating collisionless and virialized gas that makes up galactic 
halos is far from the state of gravothermal catastrophe, it is reasonable to assume that 
it is near some form of relaxation equilibrium characterized by non-extensive forms of 
entropy and energy. Admitting that the SP model follows from an idealized approach 
based on kinetic theory and an isotropic distribution function, it is interesting to verify 
empirically if the structural parameters of the halo gas can be adjusted to those of SPs, 
the equilibrium states under Tsallis' formalism. 

On the other hand, high precision N-body numerical simulations based on 
Cold Dark Matter (CDM) models, perhaps the most powerful method available for 
understanding gravitational clustering, lead to the famous results of Navarro, Frenk and 
White (NFW model) [S] that predicts density and velocity profiles which are roughly 
consistent with observations; however, it also predicts a cuspy behavior at the center of 
galaxies that is not observed in most of the rotation curves of dwarf and LSB galaxies 
O HI El EH GDI HI]- The significance of this discrepancy with observations is still under 
dispute, leading to various theoretical alternatives, either within the thermal paradigm 
(self-interacting ^2] and/or "warm" dark matter, made of lighter particles), or non- 
thermal dark matter models (with real ^21 EE GUI Id EE] or complex ^H] scalar fields, 
axions, etc), and none of these alternatives is free of controversy. However, for larger 
scales like galaxy clusters, the NFW paradigm is strongly supported by observations. A 
recent analysis 20J of ten nearby clusters using the X-ray satellite XMM-Newton has 
shown that clusters have indeed a cusped nature (see also [21] for a similar study using 
the Chandra satellite, and also [22] for a particular case study). Therefore, the CDM 
model of collision-less WIMPs remains as a viable model to account for dark matter in 
galactic halos, provided there is a mechanism to explain the discrepancies of this model 
with observations in the center of galaxies. 

Since gravity is a long-range interaction and virialized self-gravitating systems are 
characterized by non-extensive forms of entropy and energy, it is reasonable to expect 
that the final configurations of halo structure predicted by N-body simulations must be, 
somehow, related with states of relaxation associated with non-extensive formulations 
of Statistical Mechanics; therefore, a comparison between the SP and NFW models is 
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both, possible and interesting. 

A theoretical comparison between the NFW model an the SP model is established 
in [2Bj • The aim of that paper was to verify which parameters of the stellar polytropes 
provide a suitable description of the halo that resembles the one that emerges from the 
NFW model. A convenient criterion was established for a best fit of the SP model to the 
NFW profiles followed by finding the stellar polytrope whose central density, p c , central 
velocity dispersion, o~ c , and polytropic index, n, yields the same virial mass, total energy 
and maximal velocity of a given NFW halo model. Considering halo virial masses in 
the range 10 12 — 1O 15 M , it was found that the best fit to NFW profiles at all scales 
is given by polytropic indices close to n ~ 4.8, leading to an empirical estimation of 
Tsallis non-extensive parameter: q w 1.3. We will describe with more detail the main 
idea and results of that work in the present paper. 

The main purpose of our analysis is to make a dynamical analysis of two halo 
models, one based on the NFW paradigm, and other based on the SPs derived from 
Tsallis' non-extensive thermodynamics, compare them and test both with observational 
results coming from a sample of disk galaxies. 

The paper is organized as follows: in section II we provide the equilibrium 
equations of SPs and briefly summarize their connection to Tsallis' non-extensive 
entropy formalism. In section III we discuss the dynamical variables of NFW halos, 
while in section IV we describe a procedure to compare a polytropic halo with an NFW 
one and obtain numerically the parameters that characterize such polytropic halos. 
Section V deals theoretically with the dynamical consequences due to the formation 
of the galactic disk within the halo. In section VI we describe the comparison with 
observations as well as the method we used to make such a comparison. A summary of 
our results is given in section VII. 

2. Tsallis' entropy and stellar polytropes 



For a phase space given by (r, p) , the kinetic theory entropy functional associated with 
Tsallis' formalism is [21 El E] : 



where / is the distribution function and q > 1 is a real number. In the limit q — ► 1, 
the functional (0) leads to the usual Boltzmann-Gibbs functional, corresponding to 
the isothermal sphere. The condition 5 S q = leads to the distribution function that 
corresponds to the SP model characterized by the equation of state: 



where K n is a function of the polytropic index n, and can be expressed in terms of the 
central parameters: 




(1) 



p = K n p^l\ 



(2) 



K, 



n 




(3) 
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The polytropic index, n, is related to the Tsallis' parameter q > 1 by: 
3 1 

n =2 + —l (4) 
The standard approach for studying spherically symmetric hydrostatic equilibrium 
in stellar polytropes follows from inserting (j2J) into Poisson's equation, leading to the 
well known Lane-Emden equation |24| : 



with 



x 2 dx \ dx / 

/ \ Vn 



(6) 



r _ 2 47rGp c 2 p c 

1 = v r » = — • ff < = J; (7) 

G = 4.297 xlO-»^H>!, (8) 
M /kpc 

where, as mentioned above, a c and p c are the central velocity dispersion and central 
mass density respectively; G is the gravitational constant in appropriate units. Notice 
that the velocity dispersion is a measure of the kinetic temperature of the gas by means 
of the relation: a 2 = k B T c /m, with k B the Boltzmann's constant, and that we are using 
a normalization for ro, which differs from the usual one by a factor l/(n+ 1). We find 
it more convenient to consider instead of equation ((HJ) the following set of equivalent 
equilibrium equations: 

ax 

dZ n MZ l - l ' n 

dx n + 1 x 2 

where M. and Z relate to M, p (mass and mass density at radius r) by: 

z = i (10) 

Notice that in the limit n — > oo, equations become the equilibrium equations of the 
isothermal sphere. 

Once the system © has been integrated numerically, the velocity profile derived 
from the virial theorem takes the form: 

V\x) = al^ (11) 
x 

The radial distance r in kpc and enclosed mass M(r) in solar masses are given (from 
equation (|TU|)) in terms of x and M. by: 



r/kpc = 0.004220 — Q/F x, 

km/sec \ p c J 



M/M G = 944.97 f '-^-V fM^V' 2 M 
^km/sec/ \ p c J 



(12) 
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Another important dynamical quantity is the total energy of the stellar polytrope [3]: 
E = K+U 

1 4:Tir z P(r) dr - / dr — 



dr 

(13) 



leading to 



E 



1 



poly 



n — 5 



3 G M 2 ,p\i/n ,3 



(-) cr c (-(n + l)M v -{n-2)A7rr 3 Pv )],(U) 

Pc ^ 



-2 r Vc 

which must be evaluated at a fixed, but arbitrary, value of r marking a cut-off scale. 

Despite the valid general approach presented here it should be emphasized that 
there is confusion on the literature regarding the connection between n and q (equation 
(HJ)). The difference between different authors in this matter is mainly connected to 
the ambiguity on a proper definition of the one-particle distribution /(r, p) on terms 
of the probability p(r, p). In the "old" Tsallis' formalism, the relation is / oc p 
which for self-gravitating systems reveals that the extremum state for the entropy 
gives: / oc (<£>o — v 2 /2 + (^(r)) 1 /^ -1 ) where $(r) is the gravitational potential (see 
for example jS]). Such a distribution then leads to equation (J3J). The "new" Tsallis' 
formalism introduced in [2EJ establishes the "choice": / oc p q ; under this definition 
the same treatment as in leads to J2Ej: / oc ($ - v 2 /2 + ^(r))^- 1 ) and then to 
n = 1/2 + 1/(1 — q) (see also [2Z|)- Another possibility exists for the connection between 
n and q using a different approach for the treatment of self-gravitating systems [2*S] 
which leads to n = (5 — 3g)/2(l — q). There is still not a final word on which of these 
three possibilities is the most correct. 

The first two are closely related by a so called duality transformation q <-> 1/q in 
the distribution function which seems to indicate that the thermodynamics properties 
in the new formalism (/ oc p q ) can be translated into the old formalism (/ oc p) at least 
for the microcanonical ensemble description (see 26 ). 

The differences between the third description and the first two are more complicated 
and out of the scope of this article; in [22], the authors present an analytical and 
numerical study for the velocity distribution functions of self-gravitating systems and 
they claim that it is in fact this third possibility the one that better agrees with their 
results. 

Since there is yet no unique general agreement on which connection between n and 
q, and the corresponding description leading to it, is the correct one, we will use the one 
presented in equation (j3J), considering that for the purposes of this work, the different 
expressions between n and q are just different ways of parameterizing the SPs and 
therefore, using one or the other will not lead to qualitative changes in the subsequent 
results. However it is important to be clear that this might not be true and a deeper 
analysis is needed to be sure about this issue; at least, as is clear for the different 
formulas connecting n and q, the value for q found and reported in Table 1 is absolutely 
dependent on our particular choice for the description leading to equation (J3|. 
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3. NFW halos 

It is important to say that the NFW density profile [H] that we will describe in the 
present section and that will be used in this paper, is nowadays believed not to be the 
ultimate answer to the actual density profile of dark matter halos. 

As was mentioned in the Introduction, the NFW density profile shows a cuspy 
behavior in the center, in fact, for this region, it has a power law dependence that 
goes as p ~ r _1 whereas p ~ r -3 for the outer regions; such features can easily be 
seen in equation (|15|). Recent numerical simulations have shown that different profiles 
can adjust better to the simulations (see for example [SOI that propose a model where 
p ~ r -o.75 j n cen ter, see also jHIj for a recent numerical analysis and the profile 
model discussed there). Also, semianalytical models [321 CUB Ell usm g a power law 
dependence for the phase space density: p/cr 3 oc r~ a (which holds for simulated dark 
matter halos), where a is the velocity dispersion, together with the Jeans equation 
have been established a model for a better understanding of the dynamics of halos; 
predicting in particular power law dependences of the density profile going as p ~ r~ ' 8 
in the central part and p ~ r -3 44 in the external region. 

However, to simplify our analysis we will use the NFW profile as a test model due 
to its well known analytical expressions, presented in the following paragraphs, leaving 
the analysis using more accurate models as a possible future work. 

NFW numerical simulations yield the following expression for the density profile of 
virialized galactic halo structures ITTH| : 

<Wo n ^ 

1/(1 + 2/) 



where: 



Ac 3 , 

6 ° = 3[ln(l + co)-co/(l + c )]' (16) 



p = Pcrit n h 2 = 253.8 h 2 n ^, (17) 

kpc 



r 

y = c —, 



r v 

where Qq is the ratio of the total density to the critical density of the Universe, being 
1 for a flat Universe. Notice that we are using a scale parameter, y, that is different 
from that of the stellar polytropes, x. The NFW virial radius, r v , is given in terms of 
the virial mass, M v , by the condition that average halo density is proportional to the 
cosmological density p : 

A Po = —^, 19 

where A is a model-dependent numerical factor (for a ACDM model we have A ~ 100 
at z = [HE]). We remark that this last relation between the mass and virial radius 
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in terms of cosmological parameters, equation (|19|) . is valid for any halo model. The 
concentration parameter cq can be expressed in terms of M v by [3T| I38j: 

hence all quantities depend on a single free parameter M v . The mass function and 
circular velocity follow from (|15jl: 

y 



rV 3 



M = 4vr 5 oPo 
V 2 = 4nG5 oPo 



ln(l+y)- 1 + y 
Hl + y) 1 



2 



2/ 1 + 2/ 



(21) 
(22) 



The total energy of the halo (evaluated at the virial radius) is given by 

M v 2 

e nfw = - g ^t f q ( 23 ) 



F» = U(^T (24) 



where F has the approximate values: 

3 V21.5 
4. SP and NFW comparison 

In order to compare both halo models, we have to make various physically motivated 
assumptions. First, we want both models to describe a halo of the same size, but since 
the virial radius, r v , is the natural "cut-off" length scale at which the halo can be 
treated as an isolated object in equilibrium, "same size" must mean same virial mass, 
M v , by equation (JT§j) . 

Second. Both models must have the same maximum value for the rotational 
velocity obtained from equations (fTTj) and (|2~2~j) . This is a plausible assumption, as 
it is based on the Tully-Fisher relation, (35], a very well established result that has 
been tested successfully for galactic systems, showing a strong correlation between the 
total luminosity of a galaxy and its maximum rotational velocity. It can be shown 
that the Tully-Fisher relation has a cosmological origin |4TH I41j||. associated with the 
primordial power spectrum of fluctuations (the so called "cosmological Tully-Fisher 
relation"), hence it is possible to translate the correlation between maximum rotational 
velocity and total luminosity to a correlation between maximum rotational velocity and 
total (i.e. virial) mass. Since, by construction we are assuming that the polytropic and 
NFW halos have the same M v , their maximum rotational velocities must also coincide. 

Our third assumption is that the polytropic and NFW halos, complying with the 
previous requirements, also have the same total energy evaluated at the cut-off scale 
r = r v , which can easily be computed for each model. The main justification for this 
assumption follows from the fact that the total energy is a fixed quantity in the collapse 
and subsequent virialized equilibrium of dark matter halos [12]. 

|| An article is currently being prepared based on the results of |41j . 
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Since all structural variables of the NFW halo depend only on the virial mass, once 
we provide a specific value for M v all variables become determined in terms of physical 
units by means of equations (J2U| fTTH and fTH|) . Polytropic halos, on the other hand, lack 
of a closed analytic expression for mass, velocity and density profiles. In this case, 
equations (JHI) (or (JHJ)) yield numerical solutions for these profiles expressed in terms of 
the three free parameters {p c , cr c , n}. The comparison of these profiles with those of 
the NFW halos requires that we find explicit values of these free parameters, so that 
the conditions that we have outlined are met. Since we have selected three comparison 
criteria for three parameters, we have a mathematically consistent problem. 

Following the guidelines described above, we proceed to compare NFW and 
polytropic halos for M v ranging from 1O 1O M up to 10 15 M Q . From the present 
comparison we find that the values for central density, p c , of the polytropic halos are 
inversely proportional to M v , while the values for the central velocity dispersion, a c , 
are directly proportional to it (this is expected, since a c is a scale parameter in self- 
gravitating systems). The polytropic index, n, is almost constant for the selected range 
of M v , showing a very small growth as M v increases. This implies the same qualitative 
behavior of the Tsallis parameter q: it is also almost constant and is slowly increasing 
as the virial mass grows. It is worthwhile mentioning that the proportionally term K n 
in the polytropic equation of state (J2J) shows a very noticeable change, rapidly growing 
as M v increases. All these results are displayed explicitly in Table The comparison 
between both models in mass profiles, velocity profiles and density profiles is shown in 
figures HI 121 and left panel of figure El respectively. SP and NFW models have both the 
same virial mass, M v = 10 12 Mq. For other values of M v the mass, velocity and density 
profiles are qualitatively similar to the ones displayed in these figures. 

The right panel of figure El shows the slope of the logarithmic density profile as a 
function of radius. On the figure we can see more easily how the slope changes form the 
inner to the outer regions of the halo. A generalized NFW profile like the one presented 
at the beginning of section 3 (|32J |33J will lie between the SP and NFW in this 
figure. 

5. Including the galactic disk 

So far we have analyzed only the global structure of the dark matter halo without 
considering the effects of the luminous galaxy within. If one wishes to test a given 
model with observational results, it is necessary to add the baryonic disk as a dynamical 
component of the model. We present some of the general ideas of a well known method 
used to establish a relation between the dynamical properties of disk galaxies formed 
within virialized halos, with the properties of the dark matter halo itself Hoi. 

The initial conditions for the formation of the galactic disk are the following: an 
spherical halo model (SP, NFW, etc.) with a given mass, velocity and density profile; 
a total angular momentum J, total energy E and a fraction of baryionic particles 
that will constitute the disk. In a dynamical context, the formation of disks inside dark 
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matter halos can be well described with a model in which the halo responds adiabatically 
to the slow formation of the disk (due to energy dissipation of the baryonic particles). 
The method assumes that disk formation takes place in such a way that the final angular 
momentum of the disk, J d , is a fraction of the total angular momentum of the halo: 
Jd = jdJ] this transfer of angular momentum is reasonable based on a hypothesis of 
adiabatic ensemble of the disk. J d is given by: 

J d = 2tt / rV d (r)E d (r)rdr, 
Jo 

= M d h d V v / e -V-^du, (25) 
Jo V v 

where u = r/h d , V v is the virial velocity and V c (r) is the total velocity profile(disk+halo). 

The initial parameters of formation are given then by the scale quantities: M v , r v 
or V v (halo) and f d (disk), and the dynamical quantities: A = J ^ 5/2 (halo) and j d 
(disk). 

The central hypothesis of the method is that the angular momentum of each 
particle in the halo, before and after disk formation, is conserved (due to the adiabatic 
assumption). Thus, a particle that is initially at a mean radius, r i; ends up at mean 
radius, r, both related by: 

M(r;)r; = M/(r)r, (26) 

The halo evolves by several dynamical processes which end up with the baryonic 
particles collapsing into a plane to form a disk. The dark matter particles remain in 
a spherically symmetric distribution, but they do feel a gravitational pull which makes 
each mass shell of the dark halo to shrink from its original position. 

The galactic system is thus formed by two components: an spherically symmetric 
halo in virial equilibrium and a flattened disk with azimuthal symmetry on centrifugal 
equilibrium, both characterized by a radial density profile (p(r) for the halo and S(r) for 
the disk). The observed surface brightness profiles for disk galaxies favor an exponential 
form for the mass density of the baryonic disk: E(r) = H di oe~ r / hd , where h d is a scale 
radius, and E^o is the central surface density. For this distribution, the mass profile is: 

M d (r) = M d [l- e^/Mi+r/M] ? (27) 

where M d = 2iiT, d fih d . 

The rotation curve for a disk infinitely flat is given by Newtonian mechanics; a 
gravitational potential can be calculated using Poisson's equation combined with the 
hypothesis of centrifugal equilibrium and the mass distribution of the disk (equation 
(J2ZJ0; the calculation yields: 

V d 2 (r) = ^GT, dfi h d z 2 [I (z)K (z) - h^K^z)} , (28) 

where z = r/2h d . Since real galaxies are not completely flat, we will correct the last 
formula according to an observational analysis [13] > the correction diminish the value of 
V d in 5% approximately. 
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The final mass distribution of matter Mf is the sum of the dark matter inside the 
initial radius and the mass contributed by the disk: 

M f (r) = M d (r) + M(r,)(l - f d ). (29) 

The velocity profile of the halo is given by the virial theorem: 

y lu = G ( M ^; Mj(f) ) ' < 30 > 
The total rotational curve is the quadratic sum of disk and halo contributions. The 
mass distribution and the circular velocity of the disk, Md{r) and Vd(r), depend on hd 
and E^o, see (|27|28J) . The scale radius is obtained in terms of the initial conditions and 
the formula for the angular momentum of the disk (combining the expression for the 
parameter A, with (J2*5J) ): 

h d =(f)\G 1 / 2 M v r 1 J 2 \E\-^f r (31) 



where: 



(32) 



V v 

The central surface density £o,d is obtained from the Freeman relation: 

M d f d M v 

E °" " £S| " aSJ (33) 

To calculate the final mass distribution Mf we need to find the initial radius of 
contraction, r i; for each radius r, this is done by combining (J2EJ), (f2T)j) and (|2*7jl and 
solving the resulting transcendental equation. Once we obtain Mf, the velocity profile 
of the dark matter component can be calculated from equation fl3()(l . 

It is important to mention that in the end both dynamical components (disk and 
halo) depend on the value of the scale radius, h d , which turns out to depend implicitly 
on the final total circular velocity of the system (|31I32|) . but V c also depends on h d , 
so we need an iterative process to solve all the equations until a convergence for h d is 
obtained. Using the description presented, the model is complete and can be used to 
obtain total rotational curves from a set of initial conditions. 

This method can be used for different mass profiles of the halo, and allows us to 
obtain results for different halo models; in particular we will use it for the NFW and SP 
models. 

Both halo models which have been described before are characterized by isotropic 
velocity distributions for dark matter particles. However it is currently known that 
dark matter structures are actually anisotropic (see for example |3Sj and jUj, see also 
|45j . and |H)|). Such models are parametrized by the quantity j3, which is a measure 
of anisotropy: (3 = 1 — <7g/a r , where oq and a r are the tangential and radial velocity 
dispersions. In the simplest case of isotropic orbits, o~o = a r and (3 = 0. Analysis using 
numerical simulations in the past were not very clear about the value of [3: in jUj, j3 ~ 
near the center and (3 ~ 0.2 at the virial radius whereas in 0S| P ~ 0.4 at r v . More 
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recent works also suggest a universal relation between (3 and the logarithmic slope of 
the radial density distribution [33J Hlj • 

In order to include the effects of anisotropy in our analysis we would have to 
generalize the expressions for the different dynamical profiles, both for the NFW and SP 
models, using a given model for the radial dependence of (3. Also, the model described 
for the adiabatic contraction of the halo due to formation of the disk in the halo would 
have to be generalized from circular to elliptical orbits. Such analysis is out of the scope 
of the present paper. 

However, it is important to say that realistic models for dark matter halos should 
include the effects of anisotropy, then the fact that in the Tsallis' formalism described 
here, the entropy is maximized for isotropic systems tells us that such a description can 
not be the final description of the nature of dark matter structures, our intention is 
more modest, to propose it as alternative approximation to real systems. 

6. Comparison with observations 

Now that we have described a dynamical model for the galactic system (disk+halo), 
we can use it to compare the NFW and SP models against observations coming from a 
sample of real disk galaxies. 

The most direct approach would be to have observed rotation curves and surface 
brightness profiles (density distributions) for a representative sample of galaxies, then we 
could use a decomposition method to infer the components of the profile corresponding 
to dark and luminous matter for each galaxy in the sample; in this way, we would be 
able to analyze the radial distribution of each component and then, test both SP and 
NFW models to see which one fits better with the observations. However, as noticed by 
|HZj , such a sample which is also complete in galactic properties is hard to find. We are 
interested in a sample of galaxies that could be representative in the local Universe, i.e, 
that covers a wide range of luminosities, surface brightness and morphological types. 
Having such a sample will allow us to make an statistical analysis to test the SP and 
the NFW models. An alternative method to make the comparison we want is based 
in using global parameters of galaxies and the relations between them to explore the 
mass amounts of luminous and dark matter in disk galaxies. This method favors the 
statistical significance over precision. The approach used here will be on this direction, 
and we will use the outline described by [37] to make our analysis. 

In |2Zj , the authors compiled a sample with the desired characteristics: high quality 
surface photometry in the optical and near-infrared band (to obtain density profiles and 
integral luminosities), information about the rotation curves comes from the Hi or H a 
line-widths, and the total Hi gas flux to take into account the gas contribution to the 
mass distribution. This sample is composed mainly by three subsamples described in 

HOI EDI EU E2| ■ 

Some restrictions on the original samples were made in order to obtain the final 
one presented in [SZj. One of them is about the relative inclination of the galaxies: 
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the final sample of galaxies should not be highly nor slightly inclined. Otherwise either 
the surface brightness profiles in the former case, or the rotation velocities in the latter 
one, are not reliable. At the end, only galaxies with inclinations in the range of 35 to 
80 degrees were considered. Galaxies with clear signs of interaction (mainly from the 
Verjeihen sample |51j ) and with rotation curves which are still increasing at the last 
measured outer radius were excluded from the sample. Only some LSB galaxies present 
this feature. Fortunately, for most of the LSB galaxies in the sample the synthetic 
rotation curves are reported, so it was possible to perform such analysis in the exclusion 
criteria. The final sample consists of 78 galaxies and is not complete in a volume-limited 
sense, but is more or less homogeneous in the range of basic parameters: luminosity, 
surface brightness and morphological type. Dwarf galaxies were not included. 

In order to have observations from the different samples as uniform as possible, the 
raw data for each one of them was taken an then corrected for the different factors that 
alter the actual measure of a given parameter. The total magnitudes were corrected 
for galactic extinction K— correction and internal extinction jSHj- The surface 
brightness were corrected for galactic extinction, K— correction, cosmological surface 
brightness dimming, and inclination (geometrical and extinction effects). For the latter 
correction the authors in [3JJ followed the method presented in [SB] and considered 
LSB galaxies as optically thin in all bands; they defined LSB galaxies as those whose 
disk central surface brightness in the K— band after correction is larger than 18.5 
mag/arcsec 2 . The 21 cm line-widths were corrected for broadening due to turbulent 
motions and for inclination, following |51j . 

The description of the sample, homogeneous data corrections and transformations 
from observational to physical parameters is carried on by the authors in [SZ] and the 
reader is refereed to their paper for more information about it. 

With the information about the luminous disk structure parameters that we have 
from the sample (which can be used to obtain hd and £o,d)> the velocity component 
Vd(r) can be calculated usin equation (J25|) . A global quantity associated with the disk 
will be the maximum of Vd(r): 



where K m 0.193; the maximum is located at r ~ 2.2ha- The disk mass inside this 
radius is, according to equation (|2Tjl : Ma(2.2hd) = 0.64Mrf. Thus, without introducing 
further assumptions, we may define the ratio of maximum disk velocity to maximum 
total (or dynamical) velocity, Vd, m /V c , m , which is a global quantity that can be directly 
compared with theoretical predictions. 

The Vd >m /V c>m ratio is not defined at a unique radius, but it can be related to the 
dynamical, total, disk mass ratio Md yn /Md, defined at an specific radius. In particular 
at radius r m where the total rotational curve reaches its maximum we have: 




AuGKT, 04 hd 



(34) 
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where we have used the virial theorem, the Freeman relation and equation (|34|1. The 
important assumption made is that r m is proportional to the scale radius: r m = xhd- 
In fact x actually changes from galaxy to galaxy, depending on the disk surface density 
and halo mass distribution (x ~ 2.2h d for disk dominated galaxies, and x > 2.2 for halo 
dominated galaxies); is the fraction of the total disk mass at r m and depends on the 
value of x (when x = 2.2h d , f L = 0.64); however it was shown by c 37j that x changes 
slightly from galaxy to galaxy in the used sample (see figure 3a and 4a of The 
use of V dj7n /V c>m instead of M dyn /M d is suitable because it can be obtained directly from 
observational parameters, without the assumptions needed to calculate M dyn /M d . 

One of the principal results obtained in jSTj is that the ratio V± m /V c ^ m correlates 
mainly with the disk surface density Y, d . Therefore we will use this result to compare 
the NFW and the SP models with the observational results coming from the sample^f. 

To do the comparison, we need to give the values for the free parameters of the 
models. M v , r v and V v will be the same in both models, to account for same "size" 
of the halos, and they must have the same disk fraction f d , meaning that the galactic 
disk inside the spherical halos will be of the same mass; we also require the same value 
for j d , in this way the fraction of angular momentum of the disk to the halo will not 
depend, on a dynamical level, on the mass profile that the galactic system has before 
disk formation. 

We will take M v = 1 X 1O 12 M , which is a characteristic value for the mass of dark 
matter halos, as the standard value for the virial mass. For simplicity, we will assume 
that f d and j d have the same values for all disks (that is, independent of A); it turns out 
that in fact, the value of fd can not vary significantly from galaxy to galaxy, because this 
will predict larger scatter in the Tully-Fisher than the actually observed (Undergraduate 
Thesis |41j ; an article which contains these an other results is currently in preparation). 
We will also use the standard assumption in modeling disk formation: j d / f d = 1 (Ell ESI; 
although it is unclear if this hypothesis is appropriate (it depends on the efficiency of 
disk formation) we will take it as valid for the purpose of this work, arguing that the 
change of this assumption does not change the statistical trend of the presented results. 
It has not been established an appropriate value for the disk fraction f d . However, a 
plausible upper limit is the baryon fraction of the Universe as a whole (fs = fi^/fi), 
taking a Big Bang nucleosynthesis value for the baryonic density, Qb — 0.015/i~ 2 
gives an upper limit value on f d of 0.1 for a ACDM cosmological model (see also |SSj); 
however, the efficiency in forming the disk could be quite low, implying that f d could be 
substantially lower than this is confirmed by the results in [H] where an analysis on 
the scatter of the TF relation showed that f d has a numerical value even less than 0.05. 
Taking into account these considerations we will take as average values for all galaxies 

% We note that there is still controversy regarding which is the main global property of galaxies that 
correlates better with the ratio of dark to baryonic matter, some works support the results of [2Zj, see 
for example [57], and others claim that the luminosity is the key parameter, see for example |551 1551 IfiU] . 
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in both models: /<j = jj = 0.03. 

In figure |U we show a plot of the ratio Vd t m/V C:m against T, d for the described 
sample of galaxies (open circles). Despite the large scatter, an almost linear trend can 
be seen between this quantities, HSB (high surface brightness) galaxies (corresponding 
approximately to values of logS^ greater than 2.5) have greater values of Vd jm /V C)Tn than 
LSB (low surface brightness) galaxies. This means that the luminous matter content 
is greater for HSB than for LSB galaxies. The shown picture is consistent with a well 
known result: LSB galaxies are dark matter dominated systems within optical radius. 

The value of the graphic is that it allows us to bound statistically the possible 
values of the Vd tm /V Ctm ratio that galaxies with a given surface density can have. As 
was proved by [37], the size of this range of values (associated with dispersion on the 
graphic) is due to different disk mass, i.e. different virial mass, color and morphological 
type that galaxies with the same have. As has been said, the value of fa can not vary 
significantly from galaxy to galaxy, nevertheless increasing fd will increase the value of 
Vd,m/V c ^ m and viceversa, as an example, taking an interval of 0.02 around the central 
value of fa = 0.03 changes the value of the ratio Vd >m /V C)m in 4% on average, being 
less for HSB galaxies and larger for LSBs. In figure 2 of |HZj the result of varying fd 
can be seen explicitly. The observed dispersion seen in figure 0] mainly dependent on 
total disk mass and integral color (see figure 7a of [37]) has a cosmological origin which 
is explained with virial mass and concentration dispersion on the halos respectively (a 
galaxy formed inside a halo with a large value of the concentration parameter, will have 
an integral color redder than one formed inside a halo with a smaller value, an analogous 
conclusion is found for halos with different masses). 

Figure |U also shows clearly that NFW models can not reproduce at satisfaction the 
results obtained for the compiled sample without introducing unrealistic values for the 
virial mass or the concentration parameter. This is one of the results that lead us to the 
possibility of seeking an alternative to the NFW paradigm. The curves shown in figure 
H] represent both models with average values for their respective parameter; it's clear 
from the figure that the SP model follows in better agreement the average behavior of 
the observational sample than the NFW model. In terms of the parameters p c , a c and 
n of the stellar polytrope model, the scatter of the observational results is explained in 
the same way as for the NFW model, p c takes the role of the concentration parameter 
and o~ c the one corresponding to the virial mass (we found that the polytropic index n 
is almost constant for virial masses in the range 10 10 M Q to 1O 15 M ). 

In j^JJ, the same sample of observational data was also compared with results 
coming from self-consistent galaxy evolution models (see figure 3 of jHTJ)- Such semi- 
numerical models were developed in jBSl IM] ; they follow the disk galaxy formation and 
evolution in a hierarchical ACDM scenario. These models include self-consistently halo 
formation and evolution, disk star formation and feedback processes, the gravitational 
dragging of the halo due to disk formation (where the adiabatic invariance is generalized 
to elliptical orbits), secular bulge formation and other evolution processes. The overall 
main difference is a small increase on the Vd im /Vc,m ratio in figure (J3| for the more 
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realistic models compared to the simpler model described in this paper (see section 
3.2 of [37J), the difference is more notorious for HSB galaxies. It is remarkable then 
that despite the lack of real features like anisotropy, star formation and feedback, the 
simple (disk+halo) models presented above offer reasonable results compared to the 
evolutionary models, at least in what respect to figure EJ which is enough for the purpose 
of our work. 

7. Conclusions 

Motivated by the fact that stellar polytropes are the equilibrium states in Tsallis' non- 
extensive entropy formalism, we have found the structural parameters of those stellar 
polytropes that allows us to compare them with NFW halos of virial masses in the 
range 10 10 < M vir /M Q < 10 15 . The criteria for this comparison consists in demanding 
that the polytrope describes a halo having the same virial mass, virial radius, maximum 
rotational velocity and total energy as the NFW halo. These three conditions are 
sufficient to determine the three structural parameters {p c , a c , n} of the polytropic 
model; the results are displayed in Table [TJ 

We emphasize that the criteria which determine these polytropes are based on 
physically motivated assumptions: the virial radius and mass are the natural parameters 
characterizing the size of a given halo, same maximum velocity follows from the Tully- 
Fisher relation, while same total energy follows from the virilization process. As shown 
in Figure ^ the mass distribution of the polytrope grows much slower than that of the 
NFW halo up to a large radius (100 kpc) containing the core and the region where 
visible matter concentrates. Hence, as shown by Figure El the velocity profile of the 
polytrope is much less steep in the same region than that of the NFW halo. These 
features are consistent with the fact that NFW profiles predict more dark matter mass 
concentration than what is actually observed in a large sample of galaxies jHSlElEZI- 
Also, as shown in the left panel of figure El the obtained polytropes have flat cores, very 
similar to the flat isothermal cores observed in LSB galaxies (as a contrast, the cuspy 
cores of NFW halos seem to be at odds with these observations [BEl EH EH] , see also 
[HS1E])- This flat density core is a nice property, which combined with reasonable mass 
and velocity profiles, qualifies these polytropes as reasonable (albeit idealized) models 
of halo structures. 

However, in spite of their nice theoretical properties (i.e. their connection to Tsallis' 
formalism) and reasonable similarity with equivalent NFW halos, the stellar polytropes 
we have examined are very idealized configurations. Thus, we are not claiming that they 
provide a realistic description of halo structures. Instead, we suggest that their described 
features and their connection with Tsallis' formalism might indicate that the latter could 
yield useful information to understand the evolution and virialization processes of dark 
matter. Although it is necessary to pursue this idea by means of more sophisticated 
methods, including the use of numerical simulations along the lines pioneered by [3], the 
simple approach we have presented has already given interesting results. For example, 
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with respect to the parameter q( we recall that it is a free parameter of the Tsallis' 
non-extensive thermodynamics and which has not been fixed for the cosmological case), 
in this work we were able to determine its behavior as a function of the virial mass, 
and turns out to be almost constant, with values around q ~ 1.2. This result could 
be used in other contexts where extended statistical mechanics is also applied [69J. 
We have also shown that the SP model is favored over the NFW model regarding the 
dark matter content in disk galaxies (within the optical radius) which is shown by the 
average behavior of the observational sample in figure EJ These results show that the 
NFW halo model can be enhanced with the use of alternative paradigms in statistical 
mechanics, which seems to solve a recurrent item which throws a shadow in such an 
excellent description as the NFW model is. 

As mentioned, the results presented in this work show that a dark matter halo made 
out of particles which satisfy a polytropic equation of state, describes the halo in a way 
that is comparable with the description obtained from the NFW numerical simulations. 
Furthermore, our description has the advantage of not having a cuspy density profile 
near the galactic center. These results do not directly imply that dark matter halos 
obey a non-extensive entropy formalism. Further tests and experiments are needed in 
order to consider that such formalism is the one describing the dynamics of actual dark 
matter halos. At the moment, this idea is only a suggestion which is reinforced by our 
analysis. 

Finally, we would like to address the possibility of a future analysis of the 
actual statistical mechanics treatment ruling the structure formation process of dark 
matter halos by using numerical simulations to obtain the equation of state obeyed by 
such systems at different times of their evolution and then determine towards which 
equilibrium state are they moving to. Thus allowing us to have a more reliable method 
to know whether or not the Tsallis' formalism is appropriate for describing such self- 
gravitating systems. 
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Table 1. Parameters characterizing the polytropes while being compared to NFW 
halos 
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Figure 1. Mass profiles for a NFW halo with M v = 10 12 M and r v = 260 kpc (solid 
curve) and compared fit stellar polytrope (dashed curve). 
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Figure 2. Velocity profiles for the same NFW halo of figure 1 (solid line) and its 
compared stellar polytrope (dashed curve). 
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Figure 3. In the left panel the density profiles for the same NFW halo of figures 1 
and 2 (solid curve) and its compared stellar polytrope (dashed curve). In the right 
panel, the value for the slope of the logarithmic density profile is explicitly shown for 
different radii. 




Figure 4. Luminous to total dark matter content vs central surface density; open 
circles correspond to observational data, the solid line represents the NFW model and 
the dashed one the SP model. Both models have M v — 10 12 M Q . 



